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The spontaneous emission rate F of a two-level atom inside a chaotic cavity fluctuates strongly 
from one point to another because of fluctuations in the local density of modes. For a cavity with 
perfectly conducting walls and an opening containing N wavechannels, the distribution of F is given 
by P(F) oc F^/ 2 - 1 (p ^ Fo)-^-^ where Fq is the free-space rate. For small N the most probable 
value of F is much smaller than the mean value Fq. 
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The modification of the rate of spontaneous emission 
in a cavity has been a subject of extensive research i-i- 
It was shown that the cavity can both enhance and in¬ 
hibit the spontaneous emission at microwave and optical 
frequencies. The effect is due to a modification by the 
environment of the local density of modes at the position 
of the radiating atom. The efforts were concentrated on 
the fabrication of cavities of prescribed regular shape, the 
atoms being kept close to nodes or antinodes of the field 
patterns of the cavity modes. 

What can be said if the shape of the cavity is not 
regular and the exact position of the atom is unknown? 
Irregular cavities have a complicated “chaotic” field pat¬ 
tern, and it becomes difficult to state whether the spon¬ 
taneous emission rate T of a particular atom is increased 
or decreased with respect to the free-space rate Fq = 
cPuj^/STTeohc^ (corresponding to an electric dipole tran¬ 
sition with moment d, frequency ujq). Nevertheless, a 
precise statement can be made about the statistical dis¬ 
tribution of F. The distribution is universal, i.e. inde¬ 
pendent of the shape or size of the cavity, provided it is 
chaotic. 

A chaotic cavity is large compared to the wavelength 
Ao = 2ttc/ujo, and has a shape such that the light is scat¬ 
tered uniformly in phase space. (In a circular or cubic 
cavity, chaotic behavior may still occur because of diffuse 
boundary scattering or due to randomly placed scatter¬ 
ing centers.) The only parameter which enters the dis¬ 
tribution of r/Fo is the strength of the coupling of the 
cavity modes to the outside world. We assume that the 
coupling is via a hole that is small compared to the size 
of the cavity and transmits a total of N wavechannels. 
(For a hole of area A, N 27rA/AQ.) Our result for the 
distribution of F takes the universal form 


P(F) cx 


pAr/2-l 

(F + Fo)^+i’ 


( 1 ) 


shown in Fig. for several values of N. The distribution 
eventually becomes narrow and Gaussian for N 1, 
while it is still broad and strongly non-Gaussian for N as 
large as 10. The mean value of F equals Fq but the most 


probable value is smaller than Fq. 
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FIG. 1. Probability distribution of the spontaneous emis¬ 
sion rate F (normalized by the free-space rate Fo), as given 
by Eq. (|^ for several values of the number N of wavechannels 
transmitted by the hole in the cavity. For N >2 the distri¬ 
bution reaches its maximum at a rate F = Vo{N — 2)/ {N -\- 4) 
that is smaller than the mean value Fq. The variance of F 
diverges for N < 2 and equals 4Fo/(A' — 2) for larger N. The 
dashed curve is the result dil) for a hole much smaller than 
a wavelength (transmittance T = 0.1). 


As a possible experimental setup, one can imagine an 
array of cavities, each containing a few excited atoms, 
or a single cavity containing many excited atoms (widely 
separated so that they decay independently). The array 
of cavities might occur naturally in a porous material. 
Let n{t) be the number of atoms that has not decayed 
by the time t. The fraction n{t)/n{Q) is the Laplace 
transform /g°“dFP(F) exp(—Ft) of the distribution (j^), 
which is a confluent hypergeometric function. A time- 
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resolved measurenient of the emitted intensity yields n{t) 
and thereby the probability distribution P{T). Fluctu¬ 
ations of the spontaneous emission rate give rise to an 
algebraic decay n{t) oc for large t, instead of the 

usual exponential decay oc exp(—Fot). 

We proceed with the derivation of Eq. (||). We assume 
that the system is in the perturbative regime ^ , so that 
the rate of spontaneous emission is given by the Fermi 
golden rule, 


r = Aim y (AlA m{ 0 \d- E\f^) 

^ U!fj. - U!0 - i'yfj./2 


( 2 ) 


Here |0) is the initial state (excited atom -|- no photons) 
and 1 /^’'^) is the final state (atom in the ground state 
+ one photon in mode /i with frequency broadening 
7 ^). The index L or R refers to left and right eigenfunc¬ 
tions of the Maxwell equations in the open cavity, which 
form a biorthogonal set of modes. The conditions for the 
validity of perturbation theory will be discussed later. 

Eq. i) can be rewritten in terms of the local density 
of modes at the position r of the atom. 


F = 


TTUJod'^ 

- p{r,ujo), 

tlCQ 


p(r, uj) 


ii„.v 

TT 


(3) 

(4) 


where is the component along d of the electric field 
in left or right mode /r. We consider an almost empty 
cavity without any dispersive or absorptive medium in¬ 
side, in which case the distinction between the total and 
radiative density of modes @ is irrelevant. 

For a statistical description we study an ensemble of 
chaotic cavities with the same volume V and small varia¬ 
tions in shape. The average density of modes (p(r, loq)) = 
po = corresponds to the average rate Fq. Our 

aim is to find the probability distribution of p. In 
Refs. |pd]-|l^ this distribution was obtained under the 
assumption that the broadening 7 ^ was the same for all 
modes and all cavities. In our problem, the broadening 
is different for each mode and each cavity, and the dis¬ 
tribution turns out to be entirely different. 

According to the universality hypothesis of chaotic sys¬ 
tems, the statistical distribution of p can be described 
by the random-matrix theory of chaotic scattering 0. 
Starting point is the expression of the N x N scattering 
matrix S in terms of an M x M real symmetric matrix 
H (representing the discretized Helmholtz operator of the 
closed cavity) and an M x N coupling matrix W, 

S{u;) = l-2TriW\uj - H+ iTrWW'<)-'^W. (5) 


The matrix H is taken from the Gaussian orthogonal en¬ 
semble of random-matrix theory, 


P{H) cx exp [-{npoV)hr /AM] . ( 6 ) 

The limit M ^ 00 is taken at the end of the calcu¬ 
lation. The coupling matrix W has elements Wmn = 
{M/poVy^'^TT-'^dmn- 

The local density of modes is obtained from a diag¬ 
onal element of the Green function G(w) = {uj — H + 
fTrWWt)-!, 

p{frn,Uj) = -{M/nV)lmGmm{i^), (J) 

where fm is the point in space associated with the index 
TO. Because of the orthogonal invariance of P{H), the 
distribution of p is independent of to. Using Eq. (^, we 
can rewrite Eq. ( 0 ) in terms of the scattering matrix, 

p = -^^itr S'^dS/dHjnm- ( 8 ) 

ZTTl/ 

This representation of the local density of modes is the 
matrix analogue of the relationship between the local 
density of electronic states and the functional derivative 
of the scattering matrix with respect to the local electro¬ 
static potential, p(r) = {i/2Tr)tr S'^SS/SV{r). 

The matrix S'^dS/dPdmm is closely related to the ma¬ 
trix 


Q = -iS^dS/duj, (9) 

known as the Wigner-Smith time-delay matrix [^ . 
Namely, in view of Eq. (|^) we have 

itrS^dS/dHmm = Q = A^ A, (10) 

where A = . Since A is an M x matrix, the 

product AAA has M — N zero eigenvalues. The remain¬ 
ing N nonzero eigenvalues are the same as the eigenval¬ 
ues of Q, which are the so-called proper delay tim es [|l7| ] 
Ti,..., r^r. Their statistical distribution is known |lq ], 

P(ti, . . . , Tat) oc |t, - Tj I JI T--3^/2-lg-’^PoV/rfc ^ 
i<j k 

( 11 ) 

For the local density of modes (||), this implies that 


P = 


M 


N 


i=i 


( 12 ) 


where Uj is the j-th element of the eigenvector of AA^ 
corresponding to the eigenvalue Tj. In the limit M —i- 00 , 
the distribution of the vector u is Gaussian, P{u) oc 
exp(—iMIup). 

Eq. (p^, together with the distribution (11) of the 
Tj’s and the Gaussian distribution of the m^-’s, com¬ 
pletely determines the distribution of p and hence of 
F. We replace the integration over the Tj’s by the in¬ 
tegration over all elements of an arbitrary real N x 
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N matrix B such that the r^-’s are eigenvalues of 
The matrix B has distribution |18| P{B) oc 
exp(—7rpoVtrBB^)(detBB^)^^~^^^''^. Using the dimen¬ 
sionless variable x = p/=T/T q and properly rescaling 
u, B, the integral for the distribution becomes 


P{x) (X du 




X det{BB^)~^S [x — \B . 


(13) 


We first compute the distribution of the vector v = B 


which is given by Eq. (13) with the delta function re¬ 
placed by 6{v — B~^u). The result is P{v) cx (1 -I- 
Due to rotational invariance of the Gaus¬ 
sian distribution for ft, the distributions of x and |up 
are the same. Hence P{x) = f dvP(v)S(x — |up) oc 




-N-1 


This is the result (|^) announced 


in the introduction and plotted in Fig. |^. It decreases 
monotonically for N < 2, and has a maximum at non¬ 
zero T for larger N. 

This calculation holds for the so-called orthogonal sym¬ 
metry class (symmetry index /3 = 1 ), relevant for optical 
systems with time-reversal symmetry. The local density 
of states for systems with broken time-reversal symmetry 
(unitary class, /3 = 2 ) or with broken spin-rotational sym¬ 
metry (symplectic class, /3 = 4) is relevant in condensed 
matter physics. We have repeated our calculations for 
/3 = 2,4 and found = P^^(x), with p(^^(a;) given 

by Eq. ( 0 ). 

So far we have assumed that the hole in the cavity fully 
transmits at least one wavechannel, so that the trans¬ 
mittance T of the hole (the ratio of the transmitted and 
incident power) is > 1. If the hole is smaller than a wave¬ 
length, then T becomes < 1. The scattering matrix S{T) 
of the cavity coupled by a hole with transmittance T < 1 
can be expressed in terms of the scattering matrix 517 = 1 , 


S{T) = 


5|T=i + vr^ 
i + 5|T=iyr^‘ 


(14) 


To find the distribution of the local density of modes, we 
start from Eq. (||) with S replaced by S{T), repeat sim¬ 
ilar steps and average over 5 |t=i = at the end. The 
result is 


P{x) = 


xT . 


dcj) 

\/2 — T 2-\/r— T cos"^ 

[l + x{2-T + 2v^r^cos ’ 


(15) 


plotted also in Fig. 0 (dashed line, for T = 0.1). It de¬ 
creases monotonically for any P < 1 . 

The variance ((T — Tg)^) diverges if iV < 2 but the 
divergency is removed when we take into account the 
condition of applicability of the Fermi golden rule 
The perturbative treatment is valid as long as the de¬ 
cay rate T of the excited atom remains smaller than the 


width of the cavity modes contributing to the decay. 
Estimating the width of the main contributing mode as 
1 /pV = To/FpoV, we get a condition T ^ (Fo/pok’)^^^. 
Therefore, any divergent contribution of the large-F tail 
should be cutoff at F ~ (Fg/poV)^/^. The weight of the 
tail is negligibly small provided (Fg/pgV)^/^ Fg, hence 
if FgpgV = d^WgV/QTr^eglic® ^ 1. To estimate this pa¬ 
rameter, we write d = zeas {as is the Bohr radius), 
wg = 27rc/Ag, V = L^. Then FgpgV « 3.21z2a|PVAg 
is close to 1 for z = 0.17, L — 0.53 mm, Ag = 530 nm. 
We can get large room for applicability of Eqs. ( 0 ). ( 01 ) 
by going to weaker (possibly magnetic) dipoles, smaller 
cavities, or larger (possibly microwave) wavelengths. 

We conclude with a comparison with previous work 
on the local density of states in chaotic cavities EM- 
That work was motivated by different physical appli¬ 
cations (Knight shift in NMR or optical absorption). 
Our application is in a sense dual to that of Ref. [ p^ , 
where complicated electronic states interact with simple 
radiation states. Instead, we have the simplest possible 
electronic system—a two level atom—and a complicated 
structure of radiation modes. In Refs. p4|-p^ it was as¬ 
sumed that the cavity was coupled to the outside via a 
tunnel barrier of large area. In this case statistical fluc¬ 
tuations in the broadening of the levels 7 ^ (from level to 
level and from cavity to cavity) can be ignored. In the 
case of a relatively small opening, considered here, fluc¬ 
tuations of the 7 ^’s are essential. The resulting distribu¬ 
tion (0) of the local density of modes turns out to be very 
simple, compared with the result of Ref. El] (involving 
a five-fold integral in the case of unbroken time-reversal 
symmetry). We obtained our result within the frame¬ 
work of random-matrix theory. It would be interesting 
to see if it can be reproduced using the supersymmetry 
technique of Refs. 00 - 
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